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We show the existence of a region in the parameter space that defines the field dynamics in a 
Fabry-Perot cylindrical cavity, where three output stable stationary states of the light are possible 
for a given localized incident field. Two of these states do not preserve the bilateral (i.e. left-right) 
symmetry of the entire system. These broken-symmetry states are the high-transmission nonlinear 
modes of the system. We also discuss how to excite these states. 



Symmetries are at the backbone of physical theories. 
In Hamiltonian systems, symmetries are related to con- 
served quantities through Noether's theorem, and they 
provide us with a powerful tool to understand Nature. 
The left-right symmetry of the fundamental laws of Na- 
ture (apart from the nonconservation of parity for the 
weak interaction) is a well-known fact In particular, 
the electromagnetic interaction possesses this symmetry 
at the fundamental level. However, there is no reason 
why an invariance of the evolution equations should be 
an invariance of the stationary states of the system 
In fact, most of the time the state of really big systems 
does not have the symmetry of the laws which govern 
it Here in a simple symmetrical dynamical system 
with a symmetric localized driving source, we have found 
an example of the dynamical breaking of a discrete fun- 
damental symmetry of the system, namely, its left- right, 
or bilateral, symmetry. 

Shortly after the proposal and demonstration of 
nondissipative optical bistability using Fabry-Perot cavi- 
ties , the effect on the bistable behaviour of the trans- 
verse field amplitude profile of the fundamental mode of 
the cavity was taken into account Nonlinearity and 
diffraction can excite several transverse modes of the cav- 
ity, and this will considerably modify the dynamics of the 
light 0. For high finesse cavities, when only one longitu- 
dinal mode of the cavity is excited, the dynamics can be 
described with a single scalar wave equation |^J^ , which 
facilitates numerical analysis and physical interpretation. 
In this paper, we restrict ourselves to this regime. Gener- 
ally speaking, a great variety of light dynamics can be ob- 
served . Concerning bistability, it has been shown 
that diffraction gives rise to transverse instability in one 
of the two branches of the plane wave bistable regime 
This raises the question if bistability can exist at 
all when important transverse effects are taken into ac- 
count. Some families of stationary solutions have been 
found ll^, but their role in the dynamics of the light 
has not been discussed. Here we consider new families 
of stationary solutions with broken symmetry which are 
relevant to the dynamics of the light that we numerically 
observed, and we discuss their stability and excitation. 

We consider a cylindrical Fabry-Perot cavity filled 
with a self-focusing Kerr nonlinear medium, n = hq + 
n2\E\'^, where n is the refractive index of the medium. 



uq is the low-power refractive index, n2 > is the non- 
linear coefficient, and E is the electric field inside the 
cavity (see Fig. 1). The cavity is driven by a gaussian 
beam Ei with beam widths in the transverse x and y 
directions given by Wx and Wy (1/e^ half width inten- 
sity). The concave mirror in the y-direction has a ra- 
dius of curvature R, and the width Wy is chosen so that 
the incident beam matches the linear Hermite-Gaussian 
mode of the cavity in the y-direction. Both mirrors have 
very high reflectivity TZ, and the incident beam is as- 
sumed to have a Rayleigh range in the x direction = 
TTWx'^riQ/ Xq much larger than the cavity length L. Aq is 
the light wavelength in vacuum. Under these conditions, 
only a single longitudinal mode of the cavity is excited 
and the electric field E has the structure E{x,y,z,t) oc 
^(x, t) exp{—y'^/wy) sm{kcz) exp{—iujt)+c.c. , where \1/ is 
the normalized complex field envelope inside the cavity, 
Lu is the frequency of the input beam Ei, t is the time 
and kc is the resonance wavenumber of the linear cavity. 
Only one transverse mode in the y dimension is selected 
out by the use of the cylindrical cavity. The field trans- 
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where ^ d is the driving field, which is proportional to the 
incident field, ^ = x/wx is the normalized x coordinate 
and r = t/t^ is the normalized time. The normalizing 
factor to is given by = ^n^ZxIc, where c is the veloc- 
ity of light in vacuum. The normalized cavity detuning 
is given by = 2Afcz^, where AA: = kc — k and the 
wavenumber k is k = riQUojc. The amplitude decay rate 
r is F = (T + aL)zx/L, where T = 1 — 7?. is the trans- 
missivity of each mirror, and a is the loss coefficient of 
the material filling the cavity. The normalized amplitude 
q of the driving field [q = ^'^((O)), the normalized cavity 
detuning 9, and the amplitude decay rate F define the pa- 
rameter space that determines the dynamics of the light 
in the Fabry-Perot etalon. 

Let us consider a laser light that excites ^^Rb atoms 
in their D2 transition near Aq = 780 nm to provide a 
resonantly-enhanced optical nonlinearity. The mirror re- 
flectivity is TZ—Q.%%b and the cavity length L = 2 mm. 
The intracavity loss is a ~ 0.1 dB/cm and the linear 
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refractive index is taken to be ng ~ 1. A beam width 
Wx — 100 /xm corresponds to F ~ 0.15, and 9 ~ 0.4 
corresponds to a cavity frequency detuning A/ ~ 240 
MHz. The Rayleigh range in the x direction is Zx — 40 
mm. The mode beam waist in the y direction [ p^ is 
Wy — y^g^XL/nQiT{l — g)"^, where g = 1 — L/R. For 
g = 0.6, we obtain Wy ~ 28 /xm, that can be achieved 
through an appropriate cyhndrical mode-matching lens. 
The maximum nonhnear refractive index change An is 
given by An = A|5'p/47rZa;A, where A = 3/2^8 is a 
mode overlap factor given by the field profile in the y 
and z directions. For the parameters considered here, 

~ 1 corresponds to An ~ 3 x 10~^, which has been 
observed in a cylindrical Fabry-Perot experiment in ru- 
bidium vapor . 

The time evolution of the field as described by 
Eq. (|l|) , results from the interplay between different pro- 
cesses: transverse effects arising from the diffraction (sec- 
ond derivative) term in the wave equation, self-focusing 
due to the nonlinearity, energy input due to the driv- 
ing field and energy loss due to the finite transmissivity 
of both mirrors. The actual values of the parameters 
{F, 0, g}, will determine the influence on light dynamics 
of every one of the above mentioned processes. For large 
beam width, we have r,0 ^ 1, so transverse effects are 
expected to be negligible in this limit. For small beam 
width, but large enough in order to be in the regime of 
validity of Eq. , the different competing processes de- 
picted above are comparable. This is the regime we are 
considering in this paper. We will see that there is a 
region in parameter space where three possible output 
states with different broken and unbroken symmetries 
can be excited, and under appropriate circumstances, all 
stable output states can be excited. In all our consider- 
ations we will restrict ourselves to the case > 0, corre- 
sponding to a red-detuning of the laser from the cavity 
resonance. Finally, we point out that the total energy in- 
side the cavity and the transmitted beam power are both 
proportional to / = \^\^ d^, and this varies at a rate 
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Thus Eq. (|^) corresponds to a non-conservative system. 

Let us first review the main results for the stability of 
plane wave solutions against transverse perturbations Q . 
In the plane wave limit, when the second derivative of the 
field in Eq. (|l|) is neglected, the curve |5'o| — ^'d, where ^'o 
are the stationary plane wave solutions, is single- valued 
for 9 < V3r, whereas for 9 > V3F, is S-shaped and can 
lead to bistability. Stationary plane-wave solutions are 
unstable against transverse perturbation when both con- 
ditions l^-oP > r and |*oP > 6'/2 are fulfilled. Since the 
upper branch of the S-shaped curve | vE'o | — ^'d always be- 
gins at l^'oP > 2 0/3, it turns out that the upper branch 
is always unstable against transverse perturbation. This 
raises the question whether or not bistability exists at all 
when one takes into account the transverse structure of 



the beams. 

For that purpose, we first look for stationary solutions 
of Eq. (^. We set d^^/dr = 0, and solve the corre- 
sponding ordinary differential equation with a Newton- 
Raphson scheme ||l^. Note that contrary to travelling 
wave configurations, where stationary solutions are re- 
lated to a nonlinear wavenumber shift [ p^ , in the cavity 
configuration this is not the case, since the presence of 
the driving field makes Eq. (|l|) no longer invariant under 
a global phase change. We do not explore all the variety 
of stationary solutions that Eq. (|^) can support. Instead, 
we restrict ourselves to the stationary solutions that will 
play a role in our discussion of bistability. We plot in 
Fig. 2 families of stationary solutions corresponding to 
F = 0.15 and 9 = 0.4. Stationary solutions plotted in Fig. 
2 can be divided into two families. The curve labeled 
S corresponds to symmetric solutions 5'o(0, such that 
^o(C) = ^o(— 0- The curve labeled B corresponds to to 
a pair of broken-symmetry solutions, ^'o(C) and ^'o(— 0- 
The curve S has three branches defined by the sign of 
the slope of the curve in the figure, and we will refer to 
them as lower, middle, and upper branches. Figure 3 
shows the amplitude profiles of two stationary solutions, 
one with broken symmetry and the other with unbroken 
symmetry. Note the breaking of the left-right symmetry 
of the output field relative to that of the drive, for the 
broken-symmetry state. 

To investigate the stability of the stationary solu- 
tions plotted in Fig. 2, we solve Eq. (|l|) with a split- 
step Fourier Transform algorithm |l£j]. We take as in- 
put some selected slightly perturbed stationary solution 
^'(^, r = 0) = *o(C) [1 + where v{£_) is a complex 

gaussian random variable for every value of the trans- 
verse coordinate with mean fi — and typical deviation 
a = 0.01, for both the real and imaginary parts. The 
perturbative noise will seed any instability present in the 
system. Since some of the stationary solutions have bro- 
ken symmetry, in order to monitor the time evolution of 
the field, in addition to the peak amplitude, we will also 
follow the centroid defined as 
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For a symmetric beam, ^(r) = 0. After extensive nu- 
merical solutions, we find stable stationary solutions in 
the lower branch of curve S, and in curve B below 
q ~ 1.48. This means that the outcome of the simu- 
lations are the corresponding stationary solutions in 
each case. Stationary solutions corresponding to all other 
parts of curves S and B have been found to be unstable. 
In some cases, the stability analysis leads to one or the 
other of the two corresponding stable stationary solu- 
tions for a given amplitude of the driving field, which 
one depending on the noise. In other cases, the output 
is a time- varying pattern that depends on the particular 
value of the peak amplitude of the driving field. In most 
cases, when the stationary solution of the lower branch 
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of curve S is not excited, the output of simulations al- 
ways show that the beam does not preserve the bilateral 
symmetry of the driving field. For values of the peak 
amplitude of the driving field 1.37 < q < 1.48, there are 
three stable stationary solutions. The one corresponding 
to the lower branch of curve S, possesses an unbroken 
symmetry, and the transmissivity, defined as the ratio 
of the output beam power to the incident beam power, 
is low (~ 1%). The pair of solutions corresponding to 
curve B, one displaced to the left, the other to the right 
of the drive beam by equal amounts, possesses broken 
symmetry. The transmissivity of this pair of solutions 
is remarkably high (~ 75%). Thus the nonlinear Fabry- 
Perot cavity can be viewed as a nonlinear transmission 
device which selects out the broken-symmetry solutions 
as the high-transmission modes of the system. 

Numerical simulations of Eq. (|l|) show that the sym- 
metric state in the lower branch of curve S can be excited 
after switching on the driving field with the appropriate 
amplitude of the driving field. In order to excite the 
broken-symmetry states, a route to bistability has to be 
properly devised. Which one of the two broken-symmetry 
states is actually excited in the numerical simulations 
shown in Figs. 4 and 5, depends on round-off noise, and 
the method of excitation. One method is to tune the am- 
plitude of the driving field in a step-like way to a value 
q > 1.98, so that a time-varying intermediate output 
state associated with curve B is excited. Then the am- 
plitude of the driving field q is reset to the value where 
the broken-symmetry state is stable. We plot an example 
of such route to tristability in Fig. 4. Note the delay on 
the onset of the spontaneous symmetry breaking from the 
time the drive is switched down. The actual excitation 
of the broken-symmetry state depends on the character- 
istics of the intermediate unstable state that is excited, 
and on the duration of the round-trip time around the 
hysteresis loop. Another way to excite the stable state 
in curve B, it is to switch on the driving field and then 
to slowly ramp it down. Figure 5 show the actual excita- 
tion of both the symmetric and broken-symmetry states 
by this method. 

We anticipate that these results, which we have ob- 
tained at the nonlinear classical field level, should be re- 
lated to the scattering resonances of the underlying quan- 
tum field theory, in view of the correspondence principle. 
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FIG. 1. Geometry of the cylindrical Fabry-Perot cavity. 
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FIG. 2. Families of stationary solutions. Curves labelled 
S correspond to symmetric solutions, and curve labelled 
B corresponds to broken-symmetry solutions. Paxameters: 
r = 0.15 and 9 = 0.4. 
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FIG. 3. Amplitude profiel of stationary solutions from the 
S and B curves. Parameters: T = 0.15, 6 = 0.4 and q = 1.45. 

Solid line: amplitude profile of the field. Dotted line: ampli- 
tude profile of the driving field. Note the breaking of left-right 
symmetry for B. 
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FIG. 4. Step-like excitation of two stable states for 
q = 1.45. (a) Peak amplitude and (b) centroid. Parame- 
ters: F = 0.15, 9 = 0.4 and q = 1.45. S corresponds to the 
excitation of the stationary solution in the lower branch of 
curve S, B corresponds to the excitation of the stationary 
solution in curve B. Solid line: field. Dotted line: amplitude 
and centroid of the driving field. Note the delay in the onset 
of spontaneous symmetry breaking. 
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FIG. 5. Ramp-like excitation of two stable states for 
q = 1.45. (a) Peak amplitude and (b) centroid. Parame- 
ters: r = 0.15, = 0.4 and q = 1.45. S corresponds to 
the excitation of the stationary solution in the lower branch 
of curve S, and B corresponds to the excitation of the sta- 
tionary solution in curve B. Solid line: field. Dotted line: 
amplitude and centroid of the driving field. 
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